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ABSTRACT 

Context. To date, more than 30 planets have been discovered around giant stars, but only one of them has been found to be orbiting 
within 0.6 AU from the host star, in direct contrast to what is observed for FGK dwarfs. This result suggests that evolved stars 
destroy/engulf close-in planets during the red giant phase. 

Aims. We are conducting a radial velocity survey of 164 bright G and K giant stars in the southern hemisphere with the aim of studying 
the effect of the host star evolution on the inner structure of planetary systems. In this paper we present the spectroscopic atmospheric 
parameters (Teff , log^, ^, [Fe/H]) and the physical properties (mass, radius, evolutionary status) of the program stars. In addition, 
rotational velocities for all of our targets were derived. 

Methods. We used high resolution and high S/N spectra to measure the equivalent widths of many Fe i and Fe ii lines, which were 
used to derive the atmospheric parameters by imposing local thermodynamic and ionization equilibrium. The effective temperatures 
and metallicities were used, along with stellar evolutionary tracks to determine the physical properties and evolutionary status of each 
star 

Results. We found that our targets are on average metal rich and they have masses between ~ 1.0 Mq and 3.5 M©. In addition, we 
found that 122 of our targets are ascending the RGB, while 42 of them are on the HB phase. 

Key words. Stars: fundamental parameters - Stars: horizontal-branch - Planet-star interactions 



1. Introduction 

The determination of fundamental parameters (mass, radius, 
temperature, age, etc) of stars hosting exoplanets is very impor- 
tant, since it allows us to study how the physical properties of the 
host stars and the orbital parameters of the planets are related, 
when compared with non-planet host stars. This information can 
be used to test different planet formation models and to study 
the dynamical evolution of planetary systems. For instance, the 
study of the chemical abundances of stars harbouring planets led 
to the discovery of the planet-metallicity correlation for main- 
sequence stars (Gonzalez 1997; Santos et al. 2001; Fischer & 
Valenti 2005 ), which has been used as an argument in favor of 
the core-accretion model (Ida & Lin 120041 AUbert et al. 120051 
Kennedy & Kenyon|2008ll. 

During the main-sequence and the subgiant phase, the physi- 
cal parameters of stars can be derived accurately from photomet- 
ric data, since different evolutionary tracks are well separated in 
the color magnitude diagram for a given metallicity. However, 
during the red giant phase, red giant branch (RGB) and horizon- 
tal branch (HB) stars with different ages, masses and metallici- 
ties occupy a similar position in the HR diagram, making the de- 
termination of their physical parameters more difficult. In order 
to partially break this degeneracy, high-resolution spectra can be 



* Based on observations collected at La Silla - Paranal Observatory 
under programs ID's 085.C-0557 and 087.C.0476. 
** Table 1 is available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr/cgi-bin/qc at ?J/A+A/| 



used to derive effective temperature (Tgff), surface gravity (\ogg), 
and iron abundance ([Fe/H]). 

We are conducting a precision radial velocity survey of 164 
bright G and K giant stars in the southern hemisphere. The main 
goal of this project is to determine the fraction of close-in plan- 
ets (orbital periods < 150 days) orbiting RGB and HB stars, and 
compare them in order to study how the evolution of the host 
star affects the inner part of planetary systems. A more detailed 
description of the project, along with the first results will be pre- 
sented in a forthcoming paper (Jones et al. 2012 in prep.). 

In this work we present the spectroscopic atmospheric pa- 
rameters (TefF , log g , microturbulent velocity and [Fe/H]), which 
are used to derive the mass, radius and evolutionary phase of the 
program stars. In addition, rotational velocities are measured for 
our targets, so can select against rapid rotators that would pre- 
clude the measurement of precise radial velocities. 

The paper is organized as follows. In sections 2 and 3 we 
describe the targets selection, the observations and data reduc- 
tion. In section 4 we present the method used to derive the at- 
mospheric parameters. In section 5, stellar evolution models are 
used to derive the mass, radii and evolutionary status of the stars 
in our sample. In section 6 we study the dependence of macrotur- 
bulence broadening with log^ , which is used to derive projected 
rotational velocities. Finally, the summary and the discussion is 
presented in section 7. 
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Fig. 1. HR diagram including all of our targets. The red open 
circles correspond to RGB stars, while the blue filled triangles 
to HB stars. Different evolutionary tracks (Salasnich et al. l2000l l 
are overplotted for stars with 1.1 Mq, 2.0 and 3.0 M© (line 
pairs from bottom to top). The solid lines correspond to models 
with [Fe/H]=0.0 and the dashed Hnes to [Fe/H]=0.32 
(A color version of this figure is available in the online journal) 

2. Targets 

We selected a total of 164 giant stars in the Southern Hemisphere 
from the Hipparcos catalog (Ferryman el al. 1997) according 
to their position in the HR diagram (0.8 < B V < 1.2; -0.5 < 
My < 4.0). We included stars brighter than V = 8 magnitudes, 
parallaxes with a precision better than 14%, and uncertainty in 
the Johnson B -V color less than 0.02 mags. We also removed 
binary systems from the sample and we add the restriction of 
including only those stars with a Hipparcos (Hp) photometric 
variability less than 0.015 mags. The Hipparcos ID, B -V color 
and V magnitudes are listed in the first three columns of Table 
[1] The improved Hipparcos parallaxes (Van Leeuwen 2007) 
are listed in column 4, and the uncertainties are given within 
brackets. We corrected the visual magnitudes using the 3D 
extinction maps of Arenou et al. (1992) in order to compute 
absolute magnitudes and luminosities (see §5). The visual 
extinction values (Av) are listed in column 5 of Table [T] 

Figure [1] shows the position of our targets in the HR diagram 
and their resultant evolutionary status (see § 5). Evolutionary 
tracks from Salasnich et al. j2000l l are overplotted for different 
stellar masses and metallicities. In summary, our sample con- 
sists of 122 RGB stars and 42 HB stars, with a range of masses 
between -I.OMq and 3.5 Mq. 



3. Observations and data reduction. 

High resolution and high S/N spectra were taken for each of the 
stars in our sample. The targets were observed using the Fiber- 
fed Extended Range Optical Spectrograph (FEROS; Kaufer et 
al. fT999) mounted on the MPG/ESO 2.2 m telescope at La Silla 



and the Echelle Spectrograph mounted on the 1 .5 m telescope at 
CTIO, which has now been replaced by CHIRON (Schwab et al. 
|2010). FEROS has aresolving power of R ~48000, an efficiency 
of ~ 20 % and provides an almost complete optical spectral cov- 
erage (~ 3500-9200A), which allows us to study many absorp- 
tion lines in the optical range (used in the chemical analisys and 
in our radial velocity computations) and also to study the emis- 
sion in the core of the Can HK lines (3933 and 3968 A), which 
are used as chromospheric activity indicators (e.g. Jenkins et al. 
12008 ' and references therein). The echelle spectrograph mounted 
on the 1.5 m telescope, can reach a maximum resolution of R 
~ 45000 with an efficiency of ~ 1 % and covers a spectral region 
between 4020 A and 7 100 A. The exposure time of the FEROS 
targets (V < 8) ranges between 60 and 480 seconds, which leads 
to a S/N ratio between 200-300 at 5500 A and ~ 80 at 3950 A. 
The CTIO targets (V < 6) were observed with exposure times 
between 180 and 300 seconds, giving rise to S/N ratios of ~ 200 
at 5500 A. 

The FEROS spectra were reduced in a standard fashion using 
the FEROS Data Reduction Software. All the calibrations (flat- 
fields, bias and lamps) were obtained during the afternoon, ac- 
cording to the standard ESO calibration plan. The reduction of 
CTIO spectra was performed in a similar way, using an IDL- 
based pipeline available for all users. As in the case of FEROS 
data, the calibrations were taken during the afternoon, before the 
nightly stellar observations. 

4. Atmospheric parameters 

We derived spectroscopic atmospheric parameters (Teu , logg, 
[Fe/H] and microturbulent velocity) using the equivalent widths 
{EWs) of a set of neutral and singly-ionized iron lines. We used 
the 2002 version of the MOOCfl code (Sneden [T973l l. which 
solves the radiative transfer equation through a multi-layer at- 
mospheric model by imposing excitation and ionization equilib- 
rium (Saha-Boltzmann equation) and using the atomic param- 
eters for each electronic transition (excitation potential (x), os- 
cillator strength (logg/) and damping constant; see §4.1). The 
atmosphere models were obtained from the Kurucz (1993) grid. 
We linearly interpolated this grid in metallicity (fixing TgfF and 
log g ), then in Teir (fixing the metallicity and log g ) and finally 
in log g to obtain the desired atmosphere model. For a detailed 
description of this method see Gray (.2005 ) . 

4.1. Line list and atomic constants 

We adopted the line list used in Sestito et al. (I2006I I. which con- 
sists of a total of 159 Fez lines and 18 Feii lines, covering the 
spectral range between 5500A and 6800A. Features bluer than 
5500A were excluded in order to discard strongly blended lines 
and to avoid complications in the continuum tracing. Lines red- 
der than 6800 A were excluded due to the presence of many tel- 
luric features in the red part of the optical spectrum. The \oggf 
and X for each transition were also taken from Sestito et al. 
(2006; see also references therein). The collisional damping con- 
stants were computed using the Unsold (1955) approximation, 
multiplied by an enhacement factor E given by: logE = ax - b, 
where fl = 0.381 ± 0.017 and b =0.88 + 0.33. ;t coiTesponds to 
the excitation potential of the transition. This factor was derived 
from several Fez features with available accurate collisional 
damping parameters, for stars with Teif ~ 5000 K (Gratton et al. 

' http://www.as.utexas.edu/ chris/moog.html 
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Fig. 2. Comparison between our derived effective temperatures 
with those derived by three independent works. The filled black 
triangles, filled red squares and filled blue circles correspond to 
Teff 's derived by Mc William (fT990] l. Luck & Heiter (2007) and 
Randich et al. (1 19991 1. respectively. The solid line is the 1 : 1 cor- 
relation. 

(A color version of this figure is available in the online journal) 



I20031 I. Finally, we removed from the line list seven Fe i lines (AA 
5521.28, 5547.00, 5852.23, 6290.55, 6400.32, 6411. 11, 6625.04 
A) and four Feii lines {AA 5136.80, 5525.14, 5534.85, 6383.72 
A) which were found to be the most deviant points in the abun- 
dance analysis. We examined all of these lines by eye, and we 
found that all of them are heavily blended leading to an inaccu- 
rate measurement of the EW, with the exception of the feature 
at 5852.23 A, where the discrepant behavior is probably due to 
an error in its atomic parameters (log^ and/or;^). 

4.2. Equivalent widths 

Since it is very time consuming to measure EWs manually, we 
used the code ARES0(Sousa et al. l2007l l. which computes them 
automatically, by applying a gaussian fit to the spectral lines. In 
order to test the realibility of this code we also computed some 
EWs manually, using the routine splot in IRAF0. We found dif- 
ferences < 10% between our EW measurements and those de- 
rived using ARES, which are mainly due to the continuum fit. 
In the abundance analysis we included only those lines with EW 
between 10 mA and 150 mA . Lines weaker than 10 mA were 
discarded because the measured equivalent width is strongly de- 
pendent on the continuum fitting. Also, lines stronger than 150 
mA were removed from the analysis, since the gaussian fitting 



^ http://www.astro.up.pt/ sousasag/ares/ 

' IRAF is distributed by by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Fig. 3. Metallicity distribution of our targets. The red dotted line 
coiTesponds to RGB stars, while the blue dashed line to HB stars. 
The sum of both samples is shown with the solid black line. 
(A color version of this figure is available in the online journal) 



profile is not always appropiate to determine the EW and also 
because these lines are more aflrected by colhsional broadening. 

4.3. Results and uncertainties 

The eff'ective temperatures, spectroscopic surface gravities, mi- 
croturbulent velocities and iron abundances for our targets are 
presented in columns 6 - 9 of Table [T] We compared our result- 
ing effective temperatures with those derived by three other in- 
dependent studies that have a few stars in common with our sam- 
ple. These are shown in Figure |2] where we found mean differ- 
ences of (ATeff (This study - McWilliam [T990] l) = 156 + 86 K for 
14 stars in common, (ATefi (This study - Randich et al. 1999jl) 
= 108 + 70K for 6 stars in common and (ATeff (This study 
- Luck & Heiter l2007] l) = 24 + 50 K for 12 stars in common. 
Considering all of the stars in common with these three stud- 
ies we obtain (ATeff (This study - literature)) = 98 + 90 K. Based 
on this result, we adopted an uncertainty of ~ lOOK in our de- 
rived Teff 's, which is consistent with the estimated uncertain- 
ties in similar studies (see e.g. Sestito et al. 120061 Hekker & 
Melendez'200T Ghezzi et al. E&TOt . 

Figure |3] shows a histogram of the metallicity distribution of 
our targets. It can be seen that most of the giant stars in our sam- 
ple are metal rich, with ~ 50% having [Fe/H] between 0.0 and 
0.2 dex. Also, there is no significant difference in the metallic- 
ity distribution between RGB and HB stars. The uncertainties in 
the metallicities were estimated from the standard deviation of 
individual Fe i lines in the abundance analysis, and are listed in 
brackets in column 9 of Table[T] These values are larger than the 
uncertainties in the mean, but are more realistic. 

Concerning the surface gravities, several studies have re- 
vealed systematic differences between spectroscopic and pho- 
tometric logo's. In the former approach, the surface gravity is 
varied in order to obtain the same abundance for Fe i and Fe ii 



3 



M. I. Jones et al.: Study of the impact of the post-MS evolution of the host star on the orbits of close-in planets' 




Fig. 4. Difference between photometric and spectroscopic 
logg 's as a function of the effective temperature (left panel) and 
the mass (right panel) of our targets. The open red circles and 
filled blue triangles represent RGB and HB stars, respectively. 
The black dotted line shows the A log g = 0.0 boundary. 
(A color version of this figure is available in the online journal) 



lines, while in the latter case the gravity of the star is derived 
by comparing its position in the HR diagram with theoretical 
evolutionary tracks. In Figure |4] we plot the difference between 
photometric (see § 5) and spectroscopic logg 's, as a function of 
the effective temperature (left panel) and the mass (right panel) 
for our program stars. It is evident that our derived spectro- 
scopic log ^ 's are systematically lower (by ~ 0. 14 dex) than pho- 
tometric gravities for RGB stars cooler than ~ 5000 K (or less 
massive than ~ 2.0 M©). A similar discrepancy between photo- 
metric and spectroscopic \ogg 's was also found by Gratton et 
al. ( 1996) and Sestito et al. (2006). However, some authors have 
found the oppsosite result, i.e., that the spectroscopic log g 's are 
systematically higher than the photometric ones (e.g. Valenti & 
Fischer[2005 , da Silva et al. 2006). This suggests that these in- 
consistencies in the derived surface gravities are due to system- 
atic errors inherent to the method used for deriving iron abun- 
dances, mainly from Fe ii lines, which strongly affect the final 
logo's (Fei lines are quite insentive to a change in \ogg). For 
instance, an overestimation/underestimation of the position of 
the continuum, will translate into higher/lower measured EWs, 
which are matched by a higher/lower metallicity in the curve 
of growth. Since the Fe ii lines are on average weaker than Fe i 
lines, they will be more affected, hence deriving lower/higher 
logo's. However, we cannot discard a priori other effects like 
departures from LTE that might be playing a role on this (see 
e.g. Gratton et al. 1999). 

Finally, Figure |5] shows the microturbulent velocities as a 
function of the photometric \ogg 's, which are apparently more 
reliable than the spectroscopic ones, as discussed above. There 
is a clear correlation between these two parameters, where the 
microturbulent velocity decreases linearly with log g . Applying 
a linear regression yields a fit of ^ = 2.55 - 0.45 log g . A similar 
trend was also found by other authors, who obtained: ^ = 2.22 - 
0. 322 log g (Gratton et al. 1996), ^ = 1.5 - 0. 13 log g (Carretta et 
al. l2004] |. ^ =2.29 - 0.35 l og g (M onaco et al.|2005j and ^ =1 .93 - 
0.254 log g (Marino et al. l2DU8l l. 



Fig. 5. Microturbulent velocity as a function of photometric 
log g . The red open circles correspond to RGB stars and the blue 
filled triangles correspond to HB stars. A linear regression yields 
a fit of ^ = 2.55-0.45 logg, with a scatter of 0.09 kms"^ 
(A color version of this figure is available in the online journal) 

5. Masses and evolutionary status 

We used the evolutionary tracks from Salasnich et al. ( I2000I I to 
derive the mass of each star, in a similar manner to that per- 
formed in Jenkins et al. (2009) and according to the method de- 
scribed the Appendix A (only available in the online version). 
We choose to use these models because they cover a wider range 
in metallicity compared to similar evolutionary tracks (for in- 
stance Girardi et al. l2002l l. 

In addition we restricted the minimum mass of the models to 
1 Mq. This is due to the fact that less massive stars spend more 
than ~10Gyr on the main-sequence, therefore nearby (d<200 
pc), low-mass stars (< 1 M©) are not expected to have reached 
the RGB phase yet. This method allows us to derive the mass 
of a giant star given its spectroscopic (Teff and [Fe/H]) and pho- 
tometric (log L) parameters when the evolutionary status of the 
star is also known. For instance, stars cooler than ~ 5000 K and 
less luminous than ~ 1.5 Lq are ascending the RGB, and there- 
fore their masses can be derived solely by comparing their posi- 
tion in the HR diagram with stellar evolution models. However, 
more luminous RGB stars occupy a similar position in the HR 
diagram as HB stars and therefore the mass and evolutionary 
status cannot be determined simply by comparing their effective 
temperature and luminosity with isomass tracks. This is the so 
called mass-age-metallicity degeneracy, which can be partially 
broken when the metallicity of the star is known (which is the 
case for this work). 

Figure |6] shows two examples where the determination of 
the metallicity of the star is not enough to derive unambiguously 
its evolutionary status. In the upper panel we plot two models 
with the same metallicity (Z=0.008) but different masses and 
evolutionary status that cover a similar region in the HR dia- 
gram. The position of HIP21743 (Z=0.008) is also shown. It 
can be seen that this star could either be ascending the RGB 
(somewhere in between A and B) or in the HB phase (between 
C and D). However, the time scale between points C and D is 
~ 5 times longer than between A and B, therefore this star is 
most likely to be a HB star. The lower panel shows a similar sit- 
uation for HIP 68333 (Z=0.009), but this time the two models 
correspond to a I.9M0 RGB star and a I.5M0 HB giant. Both 
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Fig. 6. Upper panel: Position of HIP 2 1743 on the HR diagram 
(black dot). The two closest evolutionary models from Salasnich 
et al. (I2000I I are overplotted. The red solid line corresponds to a 
2.5 Mq RGB model, while the blue dashed line to a 2.0 Mq HB 
model. Both tracks have a metallicity of Z=0.008. Lower panel: 
same as the upper plot, but this time the masses of the models are 
1 .9 Mq for an RGB star and 1 .5M for a HB giant. The position 
of HIP 68333 is also shown. 

(A color version of this figure is available in the online journal) 



models have the same metallicity (Z=0.008). As in the previous 
case, HIP 68333 could be either somewhere in between A-B or 
C-D, considering the error bars in logTeff and LogL, but in this 
case the timescale between points C and D is ~ 45 times longer 
than between A and B. Once again, based on the timescales, this 
giant star is most likely to be on the horizontal branch. Finally, 
when the position of the star in the HR diagram is very close to 
the zero age horizontal branch, the determination of the evolu- 
tionary status is even more complicated, since the star is as likely 
to be in the red giant branch as in the He-burning phase, mainly 
due to the uncertainty in logTeff. The resulting masses, radii and 
evolutionary status of our targets are listed in columns 11, 12 
and 13 of Table [1] respectively. 



6. Rotational velocities 

Projected rotational velocities (v sin / ) were computed for the 
stars in our sample, according to the following procedure. First, 
we measured the FWHM for all of the Fei lines used in the 
abundance analysis (see § 4.1 and § 4.2) between 6000 and 6100 
A and we averaged them to obtain the total FWHM of each star 
{(Ttot)- The instrumental width {cri„s,) was measured from several 
ThAr lines present in the same spectral region, which is then 
substracted from o-,o, in order to obtain the intrinsic broaden- 
ing (cr,„,) of the star by applying: cr,-, 
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ond step consists of separating the contribution of the rotation 
and macroturbulence (non-thermal velocities), which are the two 
main line broadening mechanisms in giant stars. 



Fig. 7. Intrinsic broadening as a function of the photometric sur- 
face gravity. The red open circles correspond to RGB stars and 
the blue filled triangles to HB stars. The red open stars and blue 
filled squares correspond to RGB and HB stars from the lit- 
erature (listed in Table 2). The macroturbulence broadening is 
shown by the dashed black line. 

(A color version of this figure is available in the online journal) 



Since we have a large dataset, it can be expected that many of 
our stars show projected rotational velocities close to zero (due 
to a low inclination, i.e., sin / ~ 0) and therefore their intrinsic 
broadening is mainly due to macroturbulence broadening {(Tmac)- 
These stars can be used to determine empirically cr„,„c as a func- 
tion of some intrinsic property of the star (e.g. Melo et al. 1200 It 
Jenkins et al. 201 1). 

In order to do this, we plotted cr„„ as a function of log g , which 
is shown in Figure |7] The lower envelope (dashed black line) 
corresponds to cr„,„c, therefore stars lying close to this line are 
expected to be those with null projected rotational velocity. 
It can be seen that (Tmac increases smoothly with decreasing 
\ogg , and also seems to be systematically higher for HB stars. 
We fit a second order polynomial of the form: (Tmac = 0.2223 - 
0.0548 log g + 0.0063 log g ^, valid for RGB stars, while we just 
assumed a constant value for HB stars of cTmac =0.134 A. It is 
worth mentioning that since the macroturbulence is as depth- 
dependent phenomenon (e.g. Takeda 1995), our derived cr,„„^. is 
an average of the macroturbulence broadening for different lines, 
which are formed at different depths in the stellar atmosphere. 
Using the correlations derived above, we computed the rota- 
tional broadening (cr^-o,) by applying cr,.^, = cr,„, - o-,™,.. 

Finally, in order to convert from cr,.,,, to rotational velocities, 
we used 21 calibrators with published vsin/ 's derived by the 
Fourier transform method, which are listed in Table |2] We fit a 
straight line obtaining: vsin / = 1.18 + A2.9a-r,„, where 0-^01 is 
in A and vsin / in kmi The RMS of the fit is 0.89 kms 
We note that v sin / doesn't approach zero at 0-^0, - 0, which pro- 
duces a small systematic shift in our derived projected rotational 
velocities at low cr,.o,. We applied this conversion to all of our 
targets to finally obtain projected rotational velocities, which are 
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References, (a) Gray [19821 (b) Gray [19891 (c) Carney et aL l2008l 

listed in the last column of Table [T] 

We investigated the dependence of rotation with the lumi- 
nosity and the mass of the star Figure |8] shows our derived 
vsin/ versus logL for all of our targets. In the lower panel the 
data are binned in A log L - 0.2 dex bins, which helps to re- 
move the dispersion due to random inclination angles. It is clear 
that for RGB giants, the average rotational velocity increases 
smoothly with the luminosity. Also, it can be noticed in Figure 
[Sjthat even though HB stars rotate slightly slower than RGB gi- 
ants having the same luminosity, the difference is not statistically 
significant because of the low number of HB stars in the sample 
(~ 20 per bin). We also plotted v sin / against the mass of our tar- 
gets in Figure |9j The symbols are the same as in Figure [Hj and 
this time we binned the data in steps of A mass - 0.4 bins 
(lower panel). It can be seen that for both, RGB and HB stars, the 
average rotational velocity increases with the mass of the star, 
and no significant difference between them is observed. Also, 
this plot explains the similar trend in Figure |8] since the most 
luminous stars in the sample are also the most massive ones. 

7. Conclusions 

High resolution and high S/N spectra were used to measure 
the atmospheric parameters (Tgff , logg, v, and [Fe/H]) of 164 
giant stars, which are the targets of our precise radial velocity 
program, aimed at studying the impact of the host star evolution 
on the inner structure of planetary systems. 

We compared the resulting position in the HR diagram with 
evolutionary tracks in order to derive the physical properties of 
each star (mass, radius) and its evolutionary status. We showed 
the difficulties in the determination of the age and mass of 
stars that populate the so called "clump" in the HR diagram, 
since many evolutionary tracks for different masses and ages 
converge into this region. From the masses and radii we derived 
photometric gravities, which are systematically higher than the 
photometric values. We have also shown that the microturbulent 
velocity decreases linearly with log g , as found previously in 
other works. 
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Fig. 8. Upper panel: Projected rotational velocities against lu- 
minosity for all of our targets. The red open circles correspond 
to RGB stars while the blue filled triangles to HB stars. Lower 
panel: Same as the upper panel, but this time the data are binned 
in A log L - 0.2 dex. 

(A color version of this figure is available in the online journal) 




Mass (M, 



Fig. 9. Upper panel: Projected rotational velocities against mass 
for all of our targets. The red open circles correspond to RGB 
stars while the blue filled triangles to HB stars. Lower panel: 
Same as the upper panel, but this time the data are binned in A 
log L - 0.2 dex. 

(A color version of this figure is available in the online journal) 
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Finally we computed projected rotational velocities for all 
of our targets. We found that the rotational velocities increase 
with the mass and the luminosity of the stars, and no significant 
difl'erence is observed between RGB and HB stars. 
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Appendix A: Interpolation method 

In order to determine the mass of the star we searched for the 
closest models (isomass tracks) in the Teff - log L plane by mini- 
mizing the following quantity: 



where Teff and Lhoi are the effective temperature and the bolo- 
metric luminosity of the star. The quantities carrying an up- 
per script are the corresponding values of the model. The lu- 
minosities were derived using the bolometric corrections given 
in Alonso et al. d 19991 1. and are listed in column 10 of Table [1] 
Since the evolutionary tracks are given for different metallici- 
ties, we computed a for the two set of models containing the 
metallicity of the star (Z), i.e., < Z < Z^. Then, we used a 
linear weighting to determine the total by applying: 




(A.l) 



(A.2) 



where a and /3 correspond to the statistical weights given by: 



Z-Zi 
Z2-Z1 



1 + 



Z-Z2 
Z2-Z1 



(A.3) 



We repeated this procedure for the two isomass tracks that 
yielded the lowest total values f^^^, j and j)- We finally 
obtain the mass of the star using a linear weighting of the two 
closest solutions: 



M - kiMi + k2M2 



(A.4) 



where Mi and M2 correspond to the mass of the two isomass 
tracks. The weigthing constants are given by: 
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Table 1. Stellar parameter of the program stars 
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